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In obtaining so/utions to the first-order non/inear partial differential equations 
(PDEs) for synthesizing offset dual-shaped reflectors , it is found that previously 
observed computational problems can be avoided if the integration of the PDEs 
is started from an inner projected perimeter and integrated outward rather than 
starting from an outer projected perimeter and integrating inward . This procedure , 
however, introduces a new parameter , the main reflector inner perimeter radius p 0 , 
when given a subreflector inner angle 0 o . Furthermore , a desired outer projected 
perimeter (e.g., a circle) is no longer guaranteed. “Stability” of the integration is 
maintained if some of the initial parameters are determined first from an approxi- 
mate solution to the PDEs. A one-, two-, or three-parameter optimization algorithm 
can then be used to obtain a “best” set of parameters yielding a close fit to the de- 
sired projected outer rim. Good low cross-polarization mapping functions are also 
obtained. These methods are illustrated by synthesis of a high-gain offset-shaped 
Cassegrainian antenna and a low-noise offset-shaped Gregorian antenna. 


I. Introduction 

The problem of synthesizing dual-shaped offset re- 
flector antennas has received considerable attention in re- 
cent years, not only because it is theoretically challeng- 
ing but also because it is of considerable practical interest. 
Galindo-lsrael et al. [1] presented an exact algorithm for si- 
multaneously synthesizing the shapes of the sub and main 
reflectors for an arbitrary feed pattern, aperture ampli- 


tude, and phase distributions by solving a set of nonlinear 
first-order partial differential equations (PDEs). Jervase 
et al. [2] discussed an extension of the numerical method 
presented in [1] and gave some computed results for the 
reflector shapes. It has been observed that numerical dif- 
ficulties may arise when applying these techniques, espe- 
cially in the vicinity of the center region of the reflectors. 
This occurs when the PDEs J integration is initialized on 
the outer rims of the reflectors and proceeds “inward” [1]. 
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In this article, we investigate the possibility of cir- 
cumventing these difficulties by starting the synthesis 
procedure not from the periphery, but from inner or central 
rims of the reflectors with the integration then proceeding 
“outward” to the outer rims. The success of this proce- 
dure depends significantly on the initial conditions chosen 
at the inner rims. An approximate solution to the synthe- 
sis problem provides one good choice for the initial condi- 
tions. Numerical results for symmetric feed patterns and 
aperture distributions are included to illustrate the appli- 
cation of the synthesis procedure based on a “center” start. 
A uniform aperture distribution for a high-gain dual-offset 
shaped Cassegrainian-type reflector antenna is illustrated 
first. Next, an aperture distribution that is uniform over 
most of the aperture, but is given a Gaussian taper near 
the outer rim, is shown in a Gregorian-type configuration. 
This distribution is useful for a low-noise (but high-gain) 
large ground antenna where very low spillover is desired 
past the main reflector into a “hot” ground background. 

Starting the integration in the central region of the 
reflectors and integrating outward does not guarantee a 
circular projected outer rim as integrating inward does (if 
that is desirable). It is possible to optimize one or more 
initial parameters to obtain a desired circular outer rim 
(or other desired results). This approach was investigated 
for obtaining a circular projected outer rim with very good 
results. 

II. Numerical Method 

We choose the specified feed power pattern 1(0) and 
the desired aperture power distribution V(p) to be circu- 
larly symmetric. For the antenna system shown in Fig. 1, 
the differential equations for synthesizing the dual reflector 
system are given by [1]: 


re = r e ($,<p\p,ip,r) 

(1) 

»> = r^d, <i>\p, ip, r) 

( 2 ) 

Pe - Pe(S,<i>\p,i t ,r;p lf> ,4>^) 

( 3 ) 

'Pe = 'Pe(6,<P\p,'P,r\p^,tp 4> ) 

0 ) 


Equations (1) and (2) are obtained by applying Snell’s 
Law to the subreflector while z is calculated by using the 
path length condition viz., r -|- S — z = constant. Equa- 
tions (3) and (4) are obtained from the following three 
conditions: 

(A) Conservation of energy on a pointwise basis for geo- 
metrical optics: 

VcV(p)p\poip<j> - p<t>ip$] = 1(6) sin 9 (5) 


[PM 

Vc = J 1(0) sin OdO j J V(p)pdp (6) 

(B) Total differentiability at all points (0,<£ } p, ip) 

r$<p = (7) 

(C) Snell’s Law applied to the main reflector. 


Equations (1) through (4) are solved by starting with 
some initial values for p, 0 (very small), and ip> an d then 
increasing 0 in small increments to obtain the main and 
the subreflector surfaces in a progressive manner. The pro- 
cedure is identical to that in [1] except that the integration 
is outward (increasing 0 and p). 

III. Numerical Results for the Cassegrainian- 
Type System 

Results are presented for the following system param- 
eters: a — 11, Q — 21°, 9 m = 16°, pM = 5, r,* = 7 (see 
Fig. 1). The optical path length is assumed to be 22, the 
specified or object aperture power distribution is uniform, 
and the feed pattern is taken to be cos n 0 with n — 151 
(yielding a subreflector edge illumination level at —25.9 dB 
below the center illumination). 

As mentioned above, the reflector shapes depend 
strongly on the initial values chosen for the solution of the 
differential Eqs. (1) through (4). For a given 0 Q (minimum 
0 rim value for the subreflector), a critical initial param- 
eter is the value chosen for p 0 (minimum p rim value for 
the main reflector). As an example, Fig. 2 shows two re- 
sults for the same system parameters but with different 
initial values of po. It is observed from this figure that 
pe becomes very large near the center point (0 ^ 0) for 
one of the po choices (case B). This leads to unacceptable 
reflector shapes. 

One approach to finding good initial values is to start 
with an approximate solution (see [3], for example). An 
approximate solution can be found by relaxing the total 
differentiability condition (i.e., Eq. 7), and then prescrib- 
ing the mapping [1,3] in advance as ip = = 1, p^ = 0, 

and tpe = 0. Under these conditions the mapping from 
(9 } <j)) to (p, ip) is forced to be concentric. Thus the con- 
servation of energy reduces to 

V c V(p)ppe - 1(0) sin 0 (8) 
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The above approximate solution [3] allows a value of 
p 0 to be found for a given 9 0 which then yields very good 
and stable numerical solutions by the exact method dis- 
cussed earlier. While this value of po is not critical for 
obtaining a “good” solution, it is critical if a perfectly cir- 
cular (or other shape) is required for the outer rim of the 
main reflector. Alternate methods for choosing p 0 will be 
discussed in the next section. 

Figure 3 presents the results for ipo versus 0 for two 
cases in order to illustrate the advantage of starting from 
a central rim instead of the outer edge. Note that the out- 
ward solution (solid line) exhibits a much better behavior 
than the corresponding inward solution (dashed line). Fig- 
ure 4 shows the main and subreflector shapes derived by 
using the outward approach. Note that not only are the 
surfaces well-behaved, but also that the maximum value 
of p of the main reflector in this case does not deviate sig- 
nificantly from pMy the value for maximum p y which was 
desired to be 5. Thus, the outer periphery condition, that 
the projected aperture be circular, is approximated very 
well in this solution. 

Figure 5 shows a comparison between the radiation 
patterns in the offset plane of the shaped reflector (solid 
line) and a comparable Cassegrainian antenna. The diam- 
eter of the main reflectors of the two antennas is 80 A. The 
system parameters of the shaped reflector are A = 2200, 
Q = 23°, e M = 8.5°, PM = 1000, r* = 1000 (see Fig. 1). 
The optical units of length can be converted to wave- 
lengths by dividing by 25. The optical path length is 
assumed to be 4000. The specified main reflector aper- 
ture amplitude distribution was uniform and the feed pat- 
tern was taken to be cos” 0, with n = 151. The radiation 
patterns are calculated by using geometrical optics for the 
subreflector and physical optics for the main reflector scat- 
tering. Some radiation characteristics of the two antennas 
are given in Table 1. 


IV. Obtaining a Gregorian Design by 
Optimization of Parameters 

In this section, we will discuss a special offset dual- 
shaped Gregorian design wherein the parameter po and 
other initial parameters as well are set by an optimization 
procedure. 

The synthesis software used herein can go through a 
complete integration in less than a half minute on a 386 
20-MHz processor when approximately eight Fourier terms 
are used to represent differences in the functions 

p(<t>) and (9) 


between adjacent 6 values. Since an FFT is used for rep- 
resentation of these functions, the increase in time with 
increasing number of Fourier terms is linear. Eight Fourier 
terms was sufficient for the case to be presented. 

This very rapid two-dimensional integration of the 
nonlinear PDEs reasonably allows for the possibility of 
performing the integration many times in an optimization 
algorithm to “optimize” one or more initial parameters. 

One parameter that can be optimized is the choice 
of po for a given 9q value. Ideally, a value of po can be 
chosen at the start of the optimization by solving for the 
approximate solution to the PDEs as discussed earlier. 

The question immediately arises as to what object 
function the optimization algorithm should attempt to 
minimize. One useful object is a perfectly circular pro- 
jected perimeter of the main reflector at 0 — Omax • Thus, 
the function 


/ = SQRT 


' N 

y^pAf. - pm ) 2 /n 

n-\ 


( 10 ) 


can be minimized, where 

(PM t ' i = li •) N) 

are the perimeter values obtained at the set 
(<t>i : i = 1,...,JV) or (rpi : i = 1, N) 


(10a) 


(10b) 


for any given integration of initial parameter values. The 
value pm need not be a constant unless a circular perime- 
ter is desired. The object function could be expanded to 
obtain a prescribed set of mapping functions 


p{0,4 >) . V’OM) (U) 


as closely as possible in order to minimize cross-polarized 
currents in the aperture, for example. 

In addition to the parameter po, we have chosen 
two other initial parameters for optimization in a three- 
dimensional space. A second parameter permits one de- 
gree of freedom in the choice of the initial function 
(see [1]). For the Cassegrainian design discussed earlier, 
the value of was selected as unity. The parameter A $ 
is introduced such that 


= i + A* (^o) sin ^ ( 12 ) 
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where A y, is given in degrees. The allowable field of A y, 
is bounded to prevent caustics from occurring where not 
desired. A higher order Fourier series could be used with 
more parameters available for optimization. 

A third parameter used herein was the value of Q (see 
Fig. 1 or Fig. 6). 

Figure 6 presents the profiles of a three-dimensionally 
optimized dual-shaped Gregorian reflector pair in the off- 
set plane. A uniform phase in the main reflector aperture 
was required as indicated by the rays in Fig. 6. The object 
function was that in Eq. (10). The value 9q was chosen as 
0.1 deg, while the optimum parameters were found as 

Po = 0.015636 

Ay, = 23.53905° (13) 

and 

0 = -15.0858° 

The starting values for the above optimization were 
chosen from an approximate solution and “manual” opti- 
mization through several Alterations. Table 2 presents an 
extract of the mapping functions obtained (see Eq. 11). 

It is not clear whether this set of “optimum” parame- 
ters is global or local. Achieving a rim that is more circular 
and/or closer to the desired value of pM = 3 with better 
than three significant figures may require greater control 
of the ipy > o function (Eq. 12) than one parameter, Ay,, can 
provide. Nevertheless, the mapping is approximately cir- 
cular throughout, i.e., 

p(0,4)*p(0) ( 14 ) 

indicating that a low cross-polarization exists in the aper- 
ture. Another measure of this low cross-polarization map- 
ping is the nearly point caustic illustrated between the 
Gregorian reflectors in Fig. 6. 

Figure 7 illustrates the desired and the achieved aper- 
ture power distributions. The desired uniform phase was 
achieved exactly (see [1]). The small “hole” in the main re- 
flector near p — 0 is determined by the “optimum” choice 
of po. The difference between the geometrical optics (GO) 
computed power distribution and the desired power near 
P ~ Pm — 3.0 is at least in part due to the approximate 
GO analyses used. (The GO analyses utilized a polar grid 
wherein very wide ray “tubes” occur for larger p when 


A <j> is fixed. As A<f> was decreased, it was found that the 
discrepancy grew smaller for the same set of reflectors.) 
Similarly, and for the same reason, the small discrepancy 
near p — 0 decreased as A 9 was decreased. 

A two-dimensional optimization with Q fixed at ft — 
— 15.0858 deg, the value found with the three-dimensional 
optimization, which yields essentially the same reflector 
pair as illustrated in Fig. 6. In this case, however, the 
other values chosen for the start of the optimization were 
Po = 0.0001 and Ay, = 0.0 deg. 

Thus, we did not take advantage of the values of po 
and Ay, known from an approximate solution. The final 
values of the parameters are shown in Table 3. 

It should be noted that in order to achieve the op- 
timized values shown in Table 2 from very poor starting 
values, a careful construction of constraints had to be built 
into the optimization algorithm. However, the constraints 
were general and not limited to any particular geometry. 

Table 4 shows the resultant mapping obtained from 
the two-dimensional algorithm. It is seen that the reflec- 
tors obtained in this way differ by very little from those 
obtained and illustrated in Table 2. 

One-dimensional optimizations for an optimum value 
of po , with both Ay, and Q fixed, are generally less suc- 
cessful unless a reasonably “good” value of Ay, is chosen. 
Several cases wherein the starting value of po is chosen as: 

Po(start) = 0.0001 

are illustrated in Table 5. Values of 
Ay, = 0° 

Ay, = 10 ° 

Ay, = 20° 

were selected. When either Ay, = 0 deg or Ay, = 10 deg 
is used, a correct solution is not obtained for the PDEs. 
This is seen by the values of pm obtained at 0 = 16 deg. 
However, when Ay, — 20 deg is chosen, a good solution is 
obtained. This points up the utility of using the results 
for po and Ay, from an approximate solution as starting 
values in any optimization algorithm. 

V. Numerical Results for the Gregorian- 
Type System 

The synthesized Gregorian- type low-noise/high-gain 
dual-shaped offset reflector system was analyzed by GO 
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during the actual synthesis and with an independent GO 
computer program. The subreflector was further analyzed 
by the geometric theory of diffraction (GTD) and by phys- 
ical optics (PO).The main reflector was analyzed by PO in 
all cases. Thus, the dual-reflector analysis was performed 
by GO/PO, GTD/PO, and PO/PO (indicating the sub- 
analysis/main analysis as shown). In addition, the object 
aperture distribution, Fig. 7, was integrated directly by 
PO for comparison. 

The feed pattern used in the analysis was the same 
as that used in the synthesis. Principal polarization was 
chosen as right circular. A projection of the current dis- 
tribution on the subreflector onto a plane normal to the 
axis of the feed is shown in Fig. 8. Note the taper of 
the feed power pattern to —20 dB at the subreflector edge 
at 0 = 16 deg. This sharp taper leads to essentially no 
spillover loss past the subreflector. 

The equivalent current aperture distribution on the 
main reflector is shown in Fig. 9 when the subreflector 
is analyzed by GO. (The equivalent current is the actual 
current multiplied by the Jacobian determinant transfor- 
mation of the actual main reflector surface to the aperture 
plane.) Note that the current distribution is uniform to the 
point where it then descends Gaussian to approximately 
-20 dB at the edge of the main reflector. This distribution 
is essentially the same as the object function in Fig. 7. 

The main reflector ( D = 120A) equivalent current dis- 
tribution obtained by a GTD subreflector analysis is shown 
in Fig. 10 (a PO subreflector analysis yields little observ- 
able difference). Note the shallow diffraction ripples over 
most of the aperture (where the current was uniform in 
Fig. 9). The current distribution still has a Gaussian ta- 
per as prescribed in Fig. 7. This will result in very little 
spillover loss past the main reflector edge and thus be ideal 
for a low-noise reflector system. 

The far-field patterns (diameter of main reflector = 
120A) for each of these analysis cases is shown in Fig. 11. 
Note the very low loss in gain between the aperture in- 
tegration of the object function (Fig. 7) and the various 
diffraction analyses. This illustrates that the chosen ob- 
ject pattern has achieved the goal of potentially very low 
noise due to spillover past the main reflector into a warm 
ground environment (typical of upward-looking ground- 
based antennas). This type of aperture distribution ap- 
pears achievable only by dual-shaping methods. 


VI. Conclusions 

When integrating the nonlinear PDEs for offset dual- 
shaped reflector synthesis from a prescribed outer pro- 
jected perimeter inwards, it was found earlier that com- 
putational problems do arise in the central region of the 
reflectors [1], It is found herein that these computational 
problems are largely bypassed if the integration of the 
PDEs is started from a prescribed projected inner rim 
and integrated outward . In this case, however, a circular, 
or otherwise prescribed, outer rim is not guaranteed. In 
addition, a new parameter of the PDEs is introduced, the 
radius, po , of the inner projected rim of the main reflec- 
tor when the inner value of 0 = 6 0 for the subreflector is 
prescribed. This parameter exists in addition to other ini- 
tial parameters (e.g., ft; see Figs. 1 and 6) and the initial 
function <j > ) (see [1]). 

A poorly selected inner radius parameter p 0i or other 
poorly selected initial parameters will preclude the obtain- 
ing of useful solutions to the PDEs by the method dis- 
cussed herein or in [1]. Two complementary methods for 
choosing the parameter po and other parameters are dis- 
cussed. The first method utilizes the initial values found 
from an approximate solution to the PDEs. A second 
method utilizes an optimization algorithm which searches 
for the parameters of the PDEs that minimize a prescribed 
object function. Such an object function may measure the 
deviation of the projected outer perimeter from circular 
(or some other desired shape). Alternatively, a minimum 
cross-polarization distribution in the aperture (optimum 
mapping function) may be chosen as the object function. 

An offset dual-shaped Cassegrainian-type high-gain 
antenna is synthesized to illustrate the first method de- 
scribed above. The second method of optimization is il- 
lustrated by synthesizing a Gregorian low-noise/high-gain 
offset dual-shaped reflector antenna. One-parameter [po], 
two-parameter [po and a parameter A q, describing (<£)], 
and three-parameter [/? 0 > and ft] optimizations are il- 
lustrated. Analyses of the synthesized reflectors by GO, 
GTD, and PO are shown to verify the utility of the results. 

It should be noted that the method of using an ap- 
proximate solution to the PDEs to find the initial param- 
eters is complementary to utilization of the optimization 
algorithm since it provides an excellent set of starting val- 
ues for the parameters in the optimization procedure. 
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Table 1. Comparison of the radiation characteristics of the shaped 
reflector antenna and the Cassegralnian antenna (see Fig. 5) 


Main Reflector 

Shaped 

Reflector 

Cassegrainian 

Antenna 

Edge illumination level 

-4 dB 

-16 dB 

— 3 dB beamwidth 

0.77 deg 

0.91 deg 

1st sidelobe level 

-10 dB 

-32 dB 


Table 2. Mapping function for three-dimensional optimization, values In Eq. (13) 


Values of p at: 

4> = -90° 
4> = +90° 

^ = 0° 
i// w 180° 

<t> = +90° 

t[i - 180° 

6° 

0.156359E— 01 

0.156359E— 01 

0.156359E— 01 

0.100000E+00 

0.345170E+00 

0.353902E+00 

0.364169E+00 

0.100000E+01 

0.691 174E+00 

0.699794E+00 

0.709482E+00 

0.200000E4-01 

0.101939E+01 

0.102740E+01 

0.103634E+01 

0.300000E-I-01 

0.132247E+01 

0.132951E+01 

0.133745E+01 

0.400000E+01 

0.159466E+01 

0.160053E+01 

0.160730E+01 

0.500000E+01 

0.183220E+01 

0.183685E+01 

0.184240E+01 

0.600000E+01 

0.203348E+01 

0.203695E+01 

0.204135E+01 

0.700000E+01 

0.219891E+01 

0.220133E+01 

0.220471E+01 

0.800000E+01 

0.233064E+01 

0.23321 7E+01 

0.233470E+01 

0.900000E+01 

0.243227E+01 

0.243313E+01 

0.243500E+01 

0.100000E+02 

0.251358E+01 

0.251401E+01 

0.251552E+01 

0.110000E+02 

0.258565E+01 

0.258576E+01 

0.258707E+01 

0.120000E+02 

0.265436E+01 

0.265418E+01 

0.265537E+01 

0.130000E+02 

0.272516E+01 

0.272466E+01 

0.272582E+01 

0.140000E+02 

0.280943E+01 

0.280847E+01 

0.280969E+01 

0.150000E+02 

p M = 0.301701E+01 

0.301343E+01 

0.301529E+01 

0.160000E+02 


Table 3. Initial value parameters obtained by two-dimensional 
optimization: starting point: p 0 = 0.0001, A/,= 0.0° 


Parameters from Two-Dimensional Optimization 

Values Eq. (13) 

Values Two- 

Dimensional Optimization 

P0 = 0.015636 

p 0 = 0.015734 

A t = 23.53905° 

A ^ = 23.30467° 



Table 4, Mapping function for two-dimensional optimization, values in Table 3 



Values of p at: 



- 

4> = -90° 
ip = +90° 

0 = 0° 
iP& 0° 

<p =4-90° 
iP = 270° 

6 0 

0.157339E— 01 

0.157339E— 01 

0.157339E— 01 

O.lOOOOOE-fOO 


0.345189E+00 

0.353831E+00 

0.364174E+00 

0.100000E+01 


0.691 235E+00 

0.699701E+00 

0.709547E+00 

0.200000E+01 


0.101948E+01 

0.102728E-I-01 

0.103644E+01 

0.300000E+01 


0.132259E+01 

0.132937E+01 

0.133758E+01 

0.400000E+01 


0.159479E+01 

0.160038E+01 

0.160745E+01 

0.500000E+01 


0.183235E+01 

0.183667E+01 

0.184258E+01 

0.600000E-f 01 

= 

0.203364E+01 

0.203676E+01 

0.204154E+01 

0.700000E+01 

- 

0.219908E+01 

0.220113E+01 

0.220491 E+m 

0.800000E+01 

_ 

0.233082E-f 01 

0.233197E+01 

0.233491E+01 

0.900000E-I-01 

- 

0.243246E+01 

0.243292E+01 

0.243522E-j-01 

0.100000E+02 


0.251379E+01 

0.251378E+01 

0.251576E+01 

0.110000E+02 


0.258588E+01 

0.258550E+01 

0.258733E+01 

0.120000E+02 


0.265463E+01 

0.265387E+01 

0.265568E+01 

0.130000E+02 

_ 

0.272548E+01 

0.272429E+01 

0.272618E+01 

0.140000E+02 

- 

0.280984E+01 

0.280797E-j-01 

0.281016E+01 

0.150000E+02 


PM = 0.301 779E-J-01 

0.301343E+01 

0.301627E+01 

0.1 60000E 4-02 

= 


Table 5. Parameters and perimeter values from one-dimensional optimization of p 0 
with Ay fixed. Starting point for p 0 = 0*0001. 


<t> = -90° 
= +90° 

O 

o O 
o OO 

4> = +90° 
il> = 270° 

6° 


o 

o 

6 

II 

♦ 

Fixed 



Po = 0.019594 Final Value 


P0 = 0.195940E— 01 

0.195940E-Q1 

0.195940E— 01 

0.100000E4-00 

PM = 0.542957E+02 

0.499525E4-02 

0.459551E+02 

0.160000E4-02 


Ay = 10.0 Fixed 



Po = 0.005099 Final Value 


po = 0.509900E— 02 

0.509900E— 02 

0.509900E— 02 

0.100000E+00 

PM = 0.815349E+00 

0.697426E4-00 

0.299195E+00 

0.160000E+02 


Ay — 20.0 Fixed 



Po = 0.01747 Final Value 


po - 0.174655E— 01 

0.174655E— 01 

0.174655E-01 

0.100000E+00 

PM = 0.303909E+01 

0.300320E4-01 

0.302680E+01 

0.160000E4-02 
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Hill il 


Fig. 1. Geometry of the Cassegrainian-type dual-shaped reflector antenna system. 


120 



0, deg 



0, deg 


Fig. 3. Variation with 6 for an outward and an inward 

solution. 


Fig. 2. Variation of p and p$ with 6: (a) Initial values of po ob- 
tained from approximate solution and (b) arbitrarily chosen pQ. 





RELATIVE POWER, dB 



Fig, 4, Cassegrainian-type dual-reflector shape obtained from the synthesis procedure. 



Fig. 5. Comparison of the radiation patterns — Cassegrainlan-type 
case. 



z 


Fig. 6. Dual offset Gregorian reflectors with optimized parameters, 

POt Aip , and t C2 * 
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UNIFORM HERE 


DESIRED POWER DISTRIBUTION 
GO COMPUTED POWER DISTRIBUTION 


GAUSSIAN HERE 


Fig. 7. Main reflector aperture power distributions— desired and 
computed. 


LEGEND, CONTOUR POWER PLOT 


Fig. 8. Feed power pattern— cos 116 6 #, projected on X s 
normal to axis of feed. 
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LEGEND, CONTOUR POWER PLOT 


Fig. 9. Equivalent current aperture distribution. Subreflector anal- 
ysis by GO (see Figs. 7 and 11). 
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LEGEND, CONTOUR POWER PLOT 

Fig. 10. Equivalent current aperture distribution. Subreflector 
analysis by GTD (PQ analysis yields very similar result, see 
Fig. 11). Scale: 1 unit = 20A; radius = 60A. 


Fig. 11. Gregorlan-type dual-reflector far-field power patterns, 
where D = 120A and feed E F = cos 58 3 0 (-20 dB at edge of sub- 
reflector). Analysis by: (a) aperture PO integration (see Fig. 7); 
(b) GO for sub and PO for main (see Figs. 7 and 8); (c) GTD/PO 
(see Fig. 10); and (d) PO/PO (see Fig. 10). 





